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LARGE- AFE.TURE  SPARSE  ARRAY  ANTENNA  SYSTEMS  OF 
MODERATE  bANDWIDTH  FOR  MULTIPLE  EMITTER  LOCATION 

I.  INTRODUCTION 

The  application  of  spectral  estimation  techniques  to  resolving 
multiple,  closely-spaced  sources  of  significant  bandwidth  is  a  challenging 
problem  are?,  which  has  been  addressed  by  a  number  of  investigators  [1- 
6],  For  example,  Buckley  [1]  utilizes  linearly  constrained  minimum 
variance  beamformers  in  which  the  constraints  are  derived  from  source 
representation  spaces  and  are  termed  eigenvector  constraints.  His  paper 
contains  a  good  description  of  the  general  broadband  source  problem  plus 
analyses  of  several  approaches,  Wang  and  Kaveh  [3]  proposed  a  coherent 
subspace  approach  to  eigenvector-based  broadband  spatial  spectral  estima¬ 
tion  in  which  a  "focussed"  matrix  technique  is  employed,  based  upon 
approximate  apriori  source  location  information.  The  "focussing"  corre¬ 
sponds  to  time-delay  steering  in  broadband  beamforming.  Wax,  Shan,  and 
Kailath  [4]  utilize  decomposition  of  the  snapshot  vectors  into  a  discrete 
set  of  narrowband  frequency  components  which  are  then  processed  indepen¬ 
dently,  their  null  spectra  superimposed,  and  finally  inversion  to  form  the 
broadband  spatial  spectrum. 

This  paper  is  concerned  with  applying  a  few  of  the  above-referenced 
broadband  concepts  to  a  spatial  linear  prediction  filter  (SLPF)  source- 
estim?tor  antenna  system  which  features  a  large  aperture  width  (in 
wavelengths)  but  employs  relatively  few  aegrees-of-freedom  (DOF)  for 
accomplishing  the  estimations.  Such  a  system  would  be  a  natural  companion 

to  the  partially  adaptive,  low-sidelobe ,  phased  array  described  by  Gabriel 
Manuscript  approved  August  31,  1987. 
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[7]*  deriving  considerable  benefit  from  shared  hardware  and  data  proces¬ 
sing  equipment.  It  is  helpful  to  review  the  primary  objectives  that  led 
to  the  investigation  of  this  particular  type  of  SLPF  source  estimator 
antenna  system: 

a.  Utilize  the  largest  antenna  aperture  width  available  in  order  to 
maximize  resolution  capability  [8], 

b.  Employ  interferometer/sparse  array  techniques  to  minimize  the 
antenna  hardware  DOF  needed,  consistent  with  the  postulated  source  distri¬ 
bution  [9,10].  This  objective  is  driven  by  speed  of  response  and  cost 
considerations. 

c.  Accommodate  the  moderate  Landwidths  associated  with  phased  array 
radar  systems  [11].  This  implies  processing  received  signal  samples  which 
are  correlated  in  both  the  spatial  domain  and  the  time  domain.  Adaptive 
transversal  filters  are  often  employed  in  such  processing  [1,2,12,13],  but 
they  have  the  inherent  disadvantage  of  greatly  increasing  the  number  of 
adaptive  DOF. 

d.  Restrict  the  source  estimation  to  a  narrow  spatial  sector  that 
can  be  steered,  with  the  sector  position  determined  either  from  approxi¬ 
mate  apriori  information  or  a  deterministic  search  procedure.  This 
particular  sector  shall  be  denoted  as  the  assigned  search  (AS)  sector. 

e.  Filter  out  or  discriminate  against  all  other  sources  that  are 
located  outside  of  the  AS  sector. 

f.  Be  aminable  to  frequent  calibration  of  the  system  by  using  strong 
sources  of  opportunity. 

Section  II  develops  the  snapshot  signal  model  which  was  utilized  in 
conducting  the  investigation,  and  it  illustrates  the  problem  of  processing 


received  signal  samples  which  are  correlated  in  both  the  spatial  domain 
and  the  time  domain  over  a  modest  RF  bandwidth. 

Section  III  discusses  a  constrained  beamspace  partially  adaptive 
array  system  with  an  adaptive  transversal  filter  on  each  beam.  Here,  the 
reader  is  introduced  to  a  broadband  adaptive  system  and  becomes  acquainted 
with  dual-domain  performance  behavior  against  broadband  and  distributed 
interference  sources. 

Section  IV  then  leads  into  our  first  large-aperture,  interferometer- 
style,  SLPF  source  estimator  antenna  system  wherein  we  demonstrate  the 
estimation  of  closely  spaced  multiple  sources  of  moderate  bandwidth.  The 
SLPF  system  features  an  Interferometer  Beamformer  and  a  shaped  AS  sector 
"reference"  beam,  both  of  which  may  be  steered  to  the  AS  sector  of 
interest.  A  transversal  filter  is  connected  to  each  assigned  beam  port. 
The  initial  source  simulation  example  is  chosen  to  produce  a  wide  equiva¬ 
lent  spatial  bandwidth  spread,  but  yet  without  any  significant  spatial 
"overlap"  within  the  common  bandwidth.  This  example  is  utilized  to 
illustrate  a  number  of  principles  including;  pattern  shift  with  frequency, 
notch  cancellation  vs.  source  estimation,  and  the  very  effective  perform¬ 
ance  of  3~tap  transversal  filters  connected  into  a  few  beams. 

Section  V  considers  the  more  interesting  source  simulation  example 
wherein  the  multiple  sources  have  significant  "overlap"  of  their  equiva¬ 
lent  spatial  bandwidth  spread  and,  in  addition,  are  located  so  close 
together  as  to  require  superresolution.  This  example  is  utilized  to 
illustrate  that  simple  phase-shifter  steering  of  the  beams  may  seriously 
degrade  estimation  performance  and  that,  at  the  very  least,  the  reference 
beam  should  be  steered  via  true  time-delay.  Optimal  estimation  perform¬ 
ance  is  then  demonstrated  when  time-delay  steering  is  utilized  on  all 
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element  signals,  and  it  is  pointed  out  how  this  represents  a  practical 
Implementation  of  the  "fooussed"  covariance  matrix  technique  proposed  by 
Wang  and  Kaveh  [33 . 

It  should  be  noted  that  throughout  this  report,  source  estimation  is 
indicated  by  conventional  spatial  antenna  pattern  nulls.  This  physical 
graphical  representation  emphasizes  that  we  are  dealing  ulth  a  linear 
prediction  filter  in  which  the  null  positions  of  the  spatial  filter 
function  estimate  the  locations  of  the  spatial  point-sources.  Reference 
[8]  is  recommended  for  a  review  of  this  fundamental  principle,  and  a 
reminder  that  filter  nulls  can  often  be  rather  "fragile"  estimates  of  the 
true  spectrum. 

II.  SNAPSHOT  SIGNAL  MODEL 

This  signal  model  is  developed  in  a  tutorial  manner  because  it  is 
the  key  to  understanding  the  behavior  of  our  rather  complex,  dual-domain 
antenna  system.  Consider  a  linear  array  of  K  elements  as  shown  in  Fig.  1, 
wherein  the  received  signal  samples  are  correlated  in  both  space  and  time. 
The  postulated  signal  environment  on  any  given  observation  consists  of  I 
plane  waves  of  L  spectral  lines  arriving  from  distinct  spacial  directions 

The  RF  phase  at  the  kth  antenna  element  resulting  from  the  Jtth 
spectral  line  of  the  ith  source  will  be  the  product, 

2TTXki  Pin6i  (1) 

where  Xk^  is  the  location  of  the  element  phase  center  with  respect  to  the 
midpoint  of  the  array  in  wavelengths.  All  element  phasing  is  referenced 
to  the  midpoint  or  center  of  the  array.  Note  that  if  our  elements  are 
equally  spaced  by  a  distance,  d,  X^^  may  be  written, 
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(2> 


*»  ■  -t^)! 

where  X^  is  the  RF  wavelength  of  ith  spectral  line. 

The  complex  amplitude  of  the  1th  spectral  line  of  the  ith  source  at 
the  array  midpoint  reference  is  p^t),  such  that  we  oan  express  the  Ith 
frequency  component  of  <,ne  nth  time-^sampled  signal  at  the  kth  element  as, 

I 

Ek(n,i.)  -  nk<n,l)  ♦  l  p1(n,t)gk(ei)exp(J2xXklsine1)  (3) 

where  gk(e^)  is  the  element  pattern  amplitude  in  the  direction  and 

nk(n,l)  is  the  1th  frequency  component  of  the  nth  sample  from  the  kth 
element  independent  Gaussian  receiver  noise.  Ihe  receiver  noise  is 
assumed  to  be  a  random  prooess  with  respect  to  the  frequency  spectral  line 
index  1,  the  time  index  n  and  the  element  index  k.  It  is  also  assumed 
that  the  source  signals  are  uncorrelated  with  eceiver  noise.  Equation 
(3)  permits  us  to  construct  a  convenient  column  vector  of  observed  data  in 
the  form 

E(n,l)  -  D(l)£(n,l)  ♦  £(n,t)  (4) 

where  D(l)  is  a  Kxl  matrix  containing  a  column  vector  d^d)  for  uach  of 
the  I  source  directions,  i.e., 

dki(D  -  gk(81)exp(J2irX(<1sin0i)  .  (5) 
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Equation  (4)  separates  out  the  basic  variables  of  souroe  direction  in  the 
direction  matrix  DU),  source  baseband  signal  in  the  column  vector 
j>(n,A),  and  element  receiver  channel  noise  in  the  column  vector  £(n,A). 
The  vector  E(n,A)  is  defined  as  the  Ath  spectral  line  component  of  the 
nth  snapshot,  i.e.,  a  simultaneous  signal  sampling  across  all  K  array 
elements  at  the  nth  time  instant.  These  snapshots  would  nominally  occur 
at  the  Nyquist  sampling  rate  corresponding  to  our  receiver  bandwidth. 

The  L  frequency  spectral  line  components  of  the  source  signals  may 
be  viewed  as  a  DFT  (discrete  fourier  transform)  of  the  snapshots.  For  our 
simple  analysis  purposes,  let  us  postulate  that  the  spatial  point  source 
signals  are  uncorrelated  with  one  another  and,  in  addition,  that  the  L 
spectral  line  components  are  uncorrelated.  A  total  of  L-21  spectral  lines 
were  utilized  in  the  simulations,  equally  spaced  over  the  RF  bandwidth, 
and  the  index  A  was  stepped  sequentially  on  each  snapshot  to  produce 
"sawtooth**  swept-frsquency  signals.  Thus,  the  index  A  is  related  to  the 
index  n, 

A  -  n  -21  [lNT(|^-j]  (6) 

where  INT(x)  means  the  integer  value  of  x. 

To  illustrate  the  problem  of  processing  data  which  is  correlated  in 
both  the  frequency/time  domain  and  the  spatial  domain,  consider  forming  a 
uniform-illumination  beam  pointed  in  the  spatial  direction  0q  at  band- 
center  wavelength  For  elements  equa’ly  spaced  by  the  distance  d,  the 

element  weights,  b^,  may  be  written, 

bk  “  ^/'X:)exPi“J2lItf-^3in0o^k~^^T'^^  (7) 

o 
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The  output  of  thia  beam  for  a  single  point-source  data  veotor  constructed 
per  Eq.  (4)  would  be, 


V(n,l)-[bt.E(n,l)]-p(n,l)[bt*DU)]+[bt»rl(n,t)] 


(8) 


where  b  is  a  column  veotor  of  beam  weights  and  V(n,t)  is  the  tth  frequency 
oomponent  of  ihe  nth  snapshot  output  voltage.  The  term  of  interest  is  the 
dot  product  of  and  D(t) , 


(b^DU)] 


K 

-  I  Vk{4) 

k-1  K  K 


(9) 


Substituting  from  (5)  and  (7)  we  have  products  of  the  form 


Vk 


/  8k(8)\ 

- J«xp{jMk-*(rg4][y-«ine  -  f-sineQ]} 

\  t  to 


(10) 


The  two  domains  interact  in  the  sine  terms,  which  may  be  factored  to  yield 
the  expression, 

f 

[t- ■sine  -  T-sin'' ]-(?“)[  (Amine  -sine  ]  (11) 

t  Ao  Ao  lo 

where  fQ  is  the  band-center  RF  frequency  and  f^  is  the  tth  RF  frequency 
within  our  passband. 

The  peak  gain  of  the  uniform-illumination  beam  occurs  at  the  spatial 
direction  angle,  e,  where  Eq.  (11)  goes  to  zero,  i.e., 

f, 

( j”)sin8  -  3ineQ.  (12) 
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There  are  two  conditions  of  interest  here* 

a.  9q  -  0  deg.  The  beam  is  steered  to  boresight  and  the  direotion 
angle.  0,  will  remain  at  boresight  regardless  of  the  frequency.  Boresight 
is  defined  as  the  spatial  direction  perpendicular  to  the  linear  array 
aperture  line. 

b.  eQ*  0  deg.  The  beam  will  point  in  its  intended  steered  direc¬ 
tion,  8  .  only  when  f*  «  f  .  Otherwise,  for  f,  >  f  the  beam  shifts 

o  to  to 

toward  boresight,  and  for  f0  <  f  the  beam  shifts  away  from  boresight. 

t  o 

A  little  thought  on  Eq.  (10)  shows  that  this  shifting  vs.  frequency  occurs 
throughout  the  sidelobe  regions  of  the  beam,  i.e.,  it  affects  all  lobes 
and  nulls  except  for  the  boresight  direction.  Examples  of  this  pattern 
shifting  are  illustrated  in  the  following  sections.  The  significance  of 
off-boresight  spatial  direction  angle  shifting  vs.  frequency  is  that 
spatial  point-sources  "spread  out"  and  appear  to  the  array  as  an  equiva¬ 
lent  spatially  distributed  source  over  their  bandwidths.  This  effect  has 
serious  ramifications  in  attempting  to  estimate  the  true  spatial  positions 
of  point-sources,  and  is  the  dominant  problem  area  addressed  in  this 
paper . 

Finally,  let  us  consider  the  outputs  obtained  from  attaching  a 
tapped  delay  line  to  a  beam  output  such  as  V(n,i)  given  in  Eq.  (8).  The 
mth  tap  will  carry  the  output  voltage, 


Vra(n,A)-V(n,t)exp{-j2ir(m-1 )(-—)}  0  3) 

3 

where  f3  is  the  sampling  rate  and  Af^  is  the  translated  signal  band 
spectral  line  corresponding  to  the  4th  RF  frequency.  The  unit  delay,  T, 
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between  taps  is  assumed  to  be  exaotly  equal  to  the  sampling  periods  and  it 
is  further  assumed  that  the  signal  energy  is  band-limited  by  a  reotangular 
filter  to  prevent  "aliasing"  [14]. 

The  number  of  taps  needed  in  our  transversal  filter  may  be  estimated 
from  the  aperture-bandwidth  product  for  a  given  source  situation.  Define 
this  product  as 


B  (r-)sine,  (14) 
o 

where  D  is  the  total  antenna  aperture  width,  and  Br  is  the  RF  bandwidth 
normalized  to  the  midband  frequency.  This  product  must  be  equalled  by  the 
delay-bandwidth  product  of  the  transversal  filter,  or 

Bt(f)  -  Br(f-)ainfl1  (15) 
o 

where  t  Is  the  total  length  of  the  delay  line  and  Bfc  is  the  transversal 
filter  bandwidth  normalized  to  the  sampling  rate.  The  value  of  (t/T)  then 
gives  us  the  number  of  taps  required. 

III.  CONSTRAINED  BEAMSPACE  SLC  WITH  A  TRANSVERSAL  FILTER  ON  EACH  BEAM 

The  spatial/time/frequency  domain  adaptive  performance  is  reviewed 
for  a  constrained  beamspace  sidelobe-canceller  (SLC)  antenna  system 
wherein  each  auxiliary  SLC  beam  feeds  into  a  multiple-tap  delay  line 
(transversal  filter).  The  purpose  of  the  review  is  threefold:  (a)  to 
illustrate  the  interaction  between  the  spatial  domain  and  the  time/ 
frequency  domain  in  cancelling  broadband  spatial  point-sources;  (b)  to 
illustrate  the  special  problem  area  encountered  when  spatially  distributed 
sources  are  involved;  and  (c)  to  become  acquainted  with  an  adaptive  SLC 


9 


antenna  system  which  utilises  similar  principles  of  operation  and  would 
constitute  a  natural  application  area  for  the  postulated  source  estimation 
processor . 

Figure  2  Miustrates  a  schematic  diagram  of  the  particular  type  of 
partially  adaptive  array  system  under  discussion.  It  consists  of  two 
beamformer  subsystems  which  connect  into  common  array  aperture  elements. 
The  use  of  common  elements  helps  to  avoid  polarization  problems  in  the 
adaptive  nulling  performance.  The  mainbeam  subsystem  (on  the  left-hand 
side  in  Fig.  2)  is  intended  to  function  as  a  conventional  low-sidelobe 
phased  array  capable  of  electronic  scan  over  the  region  of  interest. 
Sidelobe  level  is  determined  by  the  usual  quality  of  the  elements,  phase 
shifter  nodules,  and  corporate  feed.  The  second  beamformer  subsystem  is 
intended  to  furnish  a  set  of  SLC  beams  for  selective  subtraction  from  the 
mainbeam.  Thus,  it  is  referred  to  as  the  SLC  subsystem.  This  SLC  sub¬ 
system  is  auxiliary  to  the  mainbeam  subsystem  and  consists  of  a  beamformer 
which  is  coupled  into  the  aperture  at  the  elements,  prior  to  the  TH  module 
phase  shifters.  This  beamformer  does  not  require  low-3idelobe  design  and 
may  consist  of  either  a  Butler  matrix  type  or  a  lens  type  such  a3  a  Rotman 
lens.  A  very  favorable  feature  from  the  standpoint  of  potential  bandwidth 
is  that  all  beams  have  the  same  phase  center,  i.e.,  the  geometric  center 
of  the  array  aperture.  Orthogonality  in  the  family  of  beams  is  desirable 
but  need  net  be  precise  for  our  purposes. 

The  output  beam  ports  connect  into  a  "Beam  Assignment  Selector" 
wherein  they  are  electronically  switched  into  the  SLC  algorithm  processor 
based  upon  source  location  estimates.  The  idea  is  to  connect  in  only 
enough  SLC  beams  to  cancel  a  given  source  distribution  situation,  thus 
minimizing  sideiobe  degradation.  The  general  principles  of  bearnspace  SLC 
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have  been  described  in  the  literature  [7]  and  shown  to  offer  advantages  of 
a  stable  raainbeam,  retention  of  low  sidelobes,  very  fast  adaptive 
response,  and  no  adaptive  grating  lobes. 

Each  of  the  selected  beams  feeds  into  a  tapped  delay  line,  usually 
referred  to  as  a  transversal  filter.  It  is  this  latter  arrangement  that 
permits  adaptivity  in  the  time/frequency  domain  [13,15,16]  in  addition  to 
adaptivity  in  the  spatial  domain.  The  transversal  filter  (TF)  output  taps 
then  feed  into  the  black  box  labelled  "SLC  Algorithm  Signal  Processor”, 
which  applies  an  adaptive  algorithm  to  obtain  tho  TF  tap  weights  for 
achieving  cancellation.  This  type  of  system  is  amenable  to  any  of  the 
current  adaptive  processing  algorithms,  including  even  analogue  versions. 

Two  different  interference  scenarios  were  selected  for  illustrating 
the  cancellation  performance  of  the  Fig.  2  SLC  antenna  system;  spacial 
point-source  cf  broad  bandwidth,  and  spatially  distributed  interference 
sources.  These  were  evaluated  via  computer  simulations  which  utilized  a 
1 6-element  linear  array  with  half -wavelength  element  spacing  as  the  common 
aperture.  The  mainbeam  subsystem  was  given  a  quiescent  Taylor  illumina¬ 
tion  taper  designed  for  30  dB  sidelobes,  and  the  SLC  subsystem  beamformer 
was  chosen  to  be  a  Butler  Matrix  type.  Figure  3  shows  a  plot  of  a  typical 
quiescent  pattern  for  this  array,  with  the  mainbeam  steered  to  -15 
degrees.  The  adaptive  weignts  for  tiie  two  scenarios  were  computed  from 
the  sample  covariance  matrix  inverse  (SMI)  algorithm  [7]. 

a.  Spacial  Point-Source  With  Bandwidth 

This  simulation  involved  a  single  33  dB  point-source  located  at 
+44  degrees  azimuth,  with  the  bandwidth  parameter  varied.  A  single  beam, 
No.  14,  from  the  Butler  matrix  SLC  beamformer  was  assigned  to  cover  this 
point-source.  Cancellation  performance  tor  this  single  beam  was  evaluated 
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with  and  without  a  transversal  filter.  Figure  4  is  a  plot  of  typical 
adapted  patterns  with  a  3~tap  transversal  filter  at  the  low,  midband,  and 
high  frequencies  of  the  10  percent  RF  bandwidth  example.  Note  the  con¬ 
siderable  shift  in  the  sidelobes  versus  frequency,  but  yet  the  null 
remains  steady  at  the  location  of  the  point-source.  Figure  5  then  shows  a 
summary  of  performance  for  this  case,  where  we  plot  the  adapted  output  in 
dB  above  receiver  noise  level  versus  the  percent  RF  bandwidth  multiplied 
by  sinie^,  where  0^  is  the  angular  location  of  the  point-source.  Note 
that  the  performance  of  the  SLC  beam  plus  a  transversal  filter  is  much 
superior  to  the  SLC  beam  alone,  which  indicates  that  the  extra  DOF's  from 
the  transversal  filter  taps  are  providing  effective  adaptivity  in  the 
frequency  domain  to  cancel  the  broadband  interference. 

The  curves  in  Fig.  5  are  dependent  upon  the  location  of  the  point- 
source  with  respect  to  both  the  mainbeam  sidelobes  and  the  auxiliary  SLC 
beam.  Thus,  they  may  be  used  only  for  qualitative  comparison  purposes. 

b.  Spatially  Distributed  Interference 

This  simulation  involved  eleven  40  dB  sources  distributed 
uniformly  in  angle  about  the  center  position  of  +44  degrees  azimuth,  with 
the  angular  width  parameter  varied.  Beams  Nos.  13*  14,  and  15  from  the 
Butler  matrix  SLC  beamformer  were  assigned  to  cover  this  distributed 
interference.  These  beams  are  illustrated  in  Fig.  6.  Cancellation  per¬ 
formance  was  evaluated  with  and  without  a  transversal  filter  on  each  beam. 
Figure  7  illustrates  a  typical  adapted  pattern  for  the  example  wherein  the 
eleven  sources  were  distributed  over  a  sector  width  equal  to  1  .25  beam- 
widths.  Note  that  the  three  beams  readily  handle  this  rather  wide  angular 
distribution  of  interference.  Figure  3  then  shows  a  summary  of  perform¬ 
ance  for  this  scenario,  where  we  plot  the  adapted  output  power  in  dB  above 
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receiver  noise  level  versus  the  distributed  sector  width  in  beamwidtrs, 
for  one,  two  and  three  SLC  beams.  Note  the  considerable  increase  in 
performance  obtained  by  utilizing  three  beams  instead  of  just  one  or  two. 

A  significant  observation  from  these  simulation  runs  was  that  the 
addition  of  a  transversal  filter  on  each  beam  had  virtually  no  effect  on 
the  cancellation  performance,  i.e.,  the  extra  DOF's  from  the  transversal 
filter  taps  did  not  provide  any  additional  adaptivity  in  the  spatial 
domain  to  counter  distributed  interference.  The  curves  in  Fig.  8  are 
dependent  upon  the  location  of  the  distributed  interference  with  respect 
to  both  the  mainbeam  sidelobes  and  the  auxiliary  SLC  beams.  Thus,  they 
may  be  used  only  for  qualitative  comparison  purposes. 

The  SLC  system  of  Fig.  2  has  also  been  found  to  be  very  effective  in 
the  cancellation  of  delayed  multipath  signals  which  get  into  the  mainbeam 
(time  domain  adaptivity)  [15]  and  the  cancellation  of  channel  mismatch 
errors  (frequency  domain  adaptivity  within  the  receiver  system  passband 
[16]. 

IV.  A  SHAPED  BEAM  SLPF  ESTIMATOR  WITH  TF 

Building  upon  the  introductory  material  of  the  previous  sections,  we 
now  proceed  to  large-aperture  interferometer /sparse  array  techniques  and 
evaluate  their  estimation  performance  against  closely-spaced  multiple 
sources  of  moderate  bandwidth.  The  particular  type  of  SLPF  source 
estimator  antenna  system  under  investigation  is  shown  in  Fig.  9.  Aperture 
width  was  chosen  to  be  25X  (wavelengths)  midband,  which  would  accommodate 
^9  elements  of  half-wavelength  spacing,  if  filled.  From  this  aperture,  we 
select  nine  (9)  elements,  spaced  3X  apart,  which  feed  their  signals  thru  a 
phase  shifter  into  an  interferometer  beamformer.  The  purpose  of  the  phase 
shifters  is  to  permit  steering  the  beamformer  beams  to  an  assigned  search 
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'.aS)  spatial  sector.  For  example,  Fig.  10  illustrates  a  typical  midband 
cluster  of  9  beams  formed  by  a  Butler  Matrix  beamformer  and  steered  to  an 
azimuth  center  position  of  42°.  The  details  of  the  interferometer  beam- 
former  weights  are  contained  in  Appendix  A,  where  it  is  noted  that  the 
cluster  is  not  unique  and  will  replicate  itself  several  times  throughout 
visible  space  because  of  the  wide  element  spacing  of  3\. 

The  output  beam  ports  connect  into  a  "Beam  Assignment  Selector” 
wherein  they  are  electronically  switched  into  a  few  transversal  filters 
(TF)  in  a  manner  similar  to  the  Fig.  2  system.  The  TF  output  taps  then 
feed  into  the  SLPF  Algorithm  Signal  Processor.  In  this  manner,  we  funnel 
down  to  the  few  DOF  that  will  be  utilized  in  performing  source  estimation 
over  a  moderate  bandwidth.  The  number  of  TF  taps  is  selectable. 

A  final  similarity  between  Fig.  9  and  Fig.  2  is  that  a  "mainbeam"  is 
included,  albeit  of  much  different  characteristics.  For  source  estima¬ 
tion,  we  desire  a  shaped  beam,  hereafter  denoted  as  the  AS  sector 
reference  beam,  which  is  of  constant  gain  across  the  AS  sector,  but  of  low 
sidelobe  level  elsewhere.  In  other  words,  the  ideal  reference  beam  would 
be  a  rectangular-shape  which  passes  the  AS  spatial  sector  but  filters  out 
all  sources  which  are  located  outside  it.  Figure  11  illustrates  the 
particular  reference  beam  that  was  utilized  for  the  simulations  contained 
in  this  report.  It  is  based  upon  a  19-element  filled  linear  array  whose 
phase  center  is  located  7.51  away  from  the  phase  center  of  the  interfero¬ 
meter  beams,  and  which  can  be  electronically  steered  to  the  AS  sector  of 
interest.  The  phase  center  of  the  reference  beam  must  be  separated 
significantly  from  the  phase  center  of  the  interferometer  beamformer  in 
order  to  achieve  robust  spatial  estimation  performance.  Details  of  the 
reference  beau  element  weights  are  contained  in  Appendix  A. 
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Equation  Formulation 


An  equation  formulation  for  the  SLPF  antenna/processor  system  of 
Fig.  9  may  be  developed  in  terms  of  the  same  pattern  subtraction 
principles  as  utilized  for  the  similar  beamspace  partially  adaptive  array 
system  of  reference  [73*  where  the  optimum  adaptive  weight  column 
vector,  W^,  was  expressed  for  K  beamformer  beams  as, 

K 

a  0  *  S*  -  I  «A  (16) 

k-1 

where  S*  is  a  column  vector  of  the  element  weights  associated  with  the  AS 
sector  reference  beam,  b  is  the  kth  Butler  Matrix  beamformer  beam  column 
vector  of  element  weights,  and  Wk  is  the  adaptive  weight  applied  to  the 
kth  beam.  The  solution  for  Wk  may  be  arrived  at  via  any  of  the  current 
adaptive  processing  algorithms,  recognizing  that  there  are  performance 
differences. 

Pointing  toward  the  sample  matrix  inverse  (SMI)  algorithm,  we  take 
advantage  of  the  fact  that  our  sample  covariance  matrix  of  signal  inputs, 

a 

R,  involves  only  the  J  assigned  beams  needed,  and  its  dimensions  reduce 
from  KxK  down  to  JxJ,  thereby  easing  the  computation  burden  involved  in 
obtaining  its  inverse  [13].  The  equivalent  "steering  vector"  A  per 
Applebaum  [17]  is  also  reduced  to  dimension  J  and  consists  of  the  cross¬ 
correlation  between  the  reference  beam  signal  V  and  J  assigned  beam 
outputs  Y, 


A  ■;  I  V(n,JOY*(n,l) 
n-1 


(17) 
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where  N  is  the  number  of  snapshots  and  index  l  is  related  to  n  (see  Eq. 
(6)).  The  jth  assigned  beam  output  for  the  nth  snapshot  signal  sample  is 
simply 


Yj (a, «. )  -  (^(n.J.)*^,),  k  set  by  j 


where  E(n,K)is  the  signal  snapshot  vector  described  in  Section  II,  and  the 
particular  beam  index  k  must  be  selected  for  the  jth  assigned  beam. 

V(n,i)  is  the  Hth  frequency  component  of  the  nth  snapshot  output  voltage 
of  the  reference  beam,  computed  from 


V(n,i)  -  (E^n.U’S*). 


vur  J  dimension  adaptive  weight  solution  thus  becomes. 


W  =  [R~\] 


where 


R-J:  I  [y(n,l)Y  t(n,t)] 
n-1 


Finally,  if  transversal  filters  with  M  taps  are  being  used,  then  the 
beam  output  snapshot  vector  becomes  a  partitioned  or  stacked  vector  of  M 
segments,  i.e., 


Y*(n,l)«[Y*(n,l)  j  Y^n.i)  j  •  ••'  Yjj(n,£)] 


wherein  the  vector  components  are  defined  as  in  Eq.  (13),  and  the  above 
expressions  must  be  modified  to  accommodate  the  partitioned  matrices. 


Needless  to  say,  ths  matrix  dimensions  are  greatly  increased  when  the 
transversal  filters  are  employed,  thus  substantially  increasing  computa¬ 
tion  burden. 

Estimating  Multiple  Sources  With  Ho  "Overlap" 

The  first  estimation  example  to  be  presented  involves  a  simulation 
in  which  we  have  three  16.7  dB  spatial  point-sources  of  10<  swept  RF 
bandwidth,  located  at  37,  42,  and  48  degrees  azimuth.  These  parameters 
were  deliberately  chosen  so  that  there  would  be  a  continuous  "equivalent” 
spatial  interference  r'.  i  <tribution  spread*  over  the  rather  wide  angular 
region  extending  from  about  35°  to  51.5°  azimuth  (roughly  5  beamwidths  in 
extent),  but  yet  without  any  significant  equivalent  spatial  overlap  within 
their  bandwldths.  To  graphically  demonstrate  the  importance  of  the  trans¬ 
versal  filters,  our  first  estimation  pattern  was  run  with  no  TF  present. 
All  nine  beams  from  the  interferometer  beamformer  were  utilized,  thus 

a 

giving  us  a  R  matrix  dimension  of  9.  Also,  both  the  AS  sector  reference 
beam  and  the  Interferometer  Beamformer  were  phased -steered  to  a  center 
position  of  42°  azimuth  as  shown  in  Figs.  10  and  11. 

A 

The  sample  covariance  matrix  R  and  the  steering  vector  A  (Eqs.  (17) 
and  (21)  were  accumulated  and  averaged  over  a  total  N-1260  snapshots.  The 
resultant  adaptive  weights  W  computed  from  Eq.  (20)  then  give  us  the 
spatial  pattern  source  estimate  shown  in  Fig.  12(a).  Several  comments  are 
in  order  regarding  this  adaptive  pattern; 

a.  There  is  no  spatial  spectrum  estimate  per  se  because,  with  no  TF 
present,  the  array  processor  cannot  separate  the  two  domains.  It  "sees"  a 


*  See  Section  II  for  an  explanation  of  equivalent  spatial  bandwidth 
spread. 


17 


distributed  source  and  utilizes  its  9  OOF  to  form  a  wide  spatial  filter 
cancellation  notch. 

b.  The  adaptive  pattern  is  shown  plotted  at  three  different  frequen¬ 
cies  within  the  10$  bandwidth)  lowest,  midband,  and  highest.  Note  that 
for  f0  >  f  the  pattern  shifts  toward  boresight,  and  for  f.  <  f  the 
pattern  shifts  away  from  boresight.  Recall  the  discussion  in  Section  II. 

c.  This  performance  is  identical  to  what  we  would  get  if  the  proces¬ 
sor  operated  from  the  9  interferometer  element  signals  directly,  i.e.,  the 
beamformer  operation  is  a  linear  transformation  of  the  element  signals. 

d.  It  demonstrates  the  same  principle  of  distributed  source  notch 
cancellation  as  shown  in  Fig.  7  for  the  SLC  system. 

A 

e.  The  eigenvalues  computed  for  the  particular  R  of  this  case, 
normalized  to  unity  receiver  noise  power  level,  were:  186.1,  160.0,  151.3, 
128.4,  112.9,  86.1,  24. 0,  2.8,  and  1.1.  Note  that  there  are  eight  unique 
eigenvalues  represented  here. 

Next,  the  same  case  was  run  with  the  TF  included,  but  with  fewer 
beams.  Beams  Nos.  3, 4, 5, 6  and  7  from  the  Interferometer  Beamformer  were 
assigned  to  cover  this  source  situation,  with  a  3“tap  TF  selected  for  each 
of  the  five  beams.  Thus,  our  matrix  dimension  for  data  processing  Js  now 
15,  and  we  must  incorporate  partitioned  vectors  per  Eq.  (22).  The 
resultant  adaptive  weights  spatial  patterns  are  shown  in  Fig.  12(b)  where 
we  note  the  following  points: 

a.  There  is  now  a  valid  spatial  spectrum  estimate  indicated  by  the 
rather  consistent  positioning  of  three  nulls  in  the  adapted  patterns, 
shown  plotted  at  three  frequencies  within  the  10$  bandwidth. 

b.  Despite  the  considerable  equivalent  bandwidth  spread,  the  3~tap 
TF  permit  separation  of  the  frequency /time  domain  from  the  spatial  domain 
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In  this  case.  Note  the  absence  of  pattern  shifting  as  compared  to  Fig. 
12(a). 

c.  This  is  an  excellent  illustration  of  a  fundamental  difference 
between  source  estimation  and  source  cancellation  systems.  Both  systems 
filter  source  signals  out  of  their  outputs,  but  in  the  former  case,  we 
desire  that  the  filtering  be  performed  via  reasonably  distinguishable 
nulls  located  at  the  source  positions  -  a  much  more  difficult  task. 

A 

d.  The  eigenvalues  computed  for  the  particular  R  of  this  case, 
normalized  to  unity  receiver  noise  power  level,  were:  *169.6,  407.5,  288.7, 
275.7,  229.1,  189.1,  157.7,  82.8,  38.5,  33.0,  9.3,  1.7,  1.6,  1.3,  and 
1.1.  Note  that  there  are  about  11  unique  eigenvalues  represented  here. 

V.  ESTIMATING  MULTIPLE  SOURCES  WITH  SIGNIFICANT  "OVERLAP* 

We  now  consider  the  more  interesting  case  of  estimating  closely 
spaced  multiple  sources  of  moderate  bandwidth  wherein  significant 
equivalent  source  overlap  exists  within  their  bandwidths.  Our  test 
simulation  example  was  formulated  by  using  the  same  three  16.7  dB  spatial 
point-sources  of  10%  swept  bandwidth  discussed  in  the  previous  section, 
but  moving  them  closer  together  to  locations  of  40,  42,  and  44  degrees 
azimuth.  Since  our  aperture  beamwidth  is  ?.1  degrees  at  these  off- 
boresight  locations,  there  is  the  added  challenge  of  superresolution, 
i.e.,  attempting  to  resolve  multiple  sources  which  are  separated  by  less 
than  a  beamwidth. 

The  SPLF  processor  sytem  conditions  are  the  same  as  for  Fig.  12(b) 
wherein  both  the  AS  sector  reference  beam  and  the  Interferometer  Beam- 
former  are  phase-steered  to  a  center  position  of  42°  azimuth  as  shown  in 
Figs.  10  and  11.  Also,  five  beams,  Nos.  3,  4,  5,  6,  and  7,  are  assigned 
to  cover  the  source  situations  with  a  3~tap  TF  selected  for  each  of  the 
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five  beans.  The  i  aultant  adaptive  weights  spatial  patterns  are  shown  in 
Fig.  13(a)  where  we  note  the  formation  of  a  good  cancellation  filter  notch 
for  the  three  sources,  but.  no  valid  spatial  spectrum  estimate. 

The  failure  to  resolve  the  sources  is  caused  by  two  factors  whioh 
were  simultaneously  present:  (a)  equivalent  spatial  domain  overlap  within 
the  common  bandwidth,  and  (b)  the  AS  sector  reference  beam  undergoes  both 
pattern  shifting  and  phase-center  shifting  vs.  frequency.  The  estimator 
system  with  TF  can  tolerate  one  of  these  factors  at  a  time,  as  in  Fig. 
12(b),  but  not  both.  To  demonstrate  this  point  with  the  second  option, 
let  us  next  remove  the  equivalent  source  shifting  in  the  AS  sector 
reference  beam  by  employing  true  time-delay  steering  instead  of  phase¬ 
steering. 

Time-delay  steering  produces  element  weights  which  vary  with 
frequency.  hecalling  the  discussion  in  Section  II,  it  is  readily  seen 
that  Eq.  (7)  would  become 

bkU)  -  (|/:^:)exp{-j2Tr(|-)3in0olk-(^-)]}  (23) 

N 

where  we  note  the  frequency  dependence  thru  1^.  Substituting  this  change 
into  Eq.  (10),  we  aee  that  there  would  be  no  beam  shifting  vs.  frequency 
in  the  steering  direction  9q.  Furthermore,  the  shifting  will  remain 
minimal  for  directions  within  a  beamwidth  of  6q,  provided  the  RF  bandwidth 
is  modest.  Appendix  A  contains  pattern  plots  of  the  AS  sector  reference 
beam  for  both  phase-steering  and  true  time-delay  steering  to  a  center 
position  of  42°  azimuth. 

This  "fixing"  or  stabilization  of  the  multiple  swept  sources  in  the 
reference  beam  can  be  better  appreciated  via  close  examination  of  Eq.  (17) 
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in  the  previous  Section  IV,  where  we  see  th;*t  the  equivalent  steering 
veotor,  A,  oonsists  of  the  cross-oor relations  between  the  reference  beam 
signals,  V(n,l),  and  the  inter feror.' ter  beam  outputs  Y(n,l).  Obviously, 
if  the  sources  remain  constant  in  the  V(n,l)  signal  with  respect  to  the 
interferometer  array  phase  center,  then  all  of  the  amplitude /phase  varia¬ 
tions  vs.  frequency  will  be  confined  to  the  interferometer  beam  outputs 
and  oan  be  compensated  for  by  the  adaptive  transversal  filters.  Figure 
13(b)  demonstrates  the  effective  compensation  by  the  TF  under  this 
condition  wherein  the  reference  beam  is  steered  via  trre  time-delay  rather 
than  via  phase  shifters.  Note  that  we  now  have  a  valid  superresolution 
spatial  spectrum  estimate  of  the  multiple  swept-bandwidth,  "overlapped" 
sources  and,  furthermore,  that  there  is  very  little  pattern  shifting 

throughout  the  entire  AS  sector  (compare  against  Fig.  12(b)). 

* 

The  eigenvalues  computed  for  the  particular  R  of  this  case, 
normalized  to  unity  receiver  noise  power  level,  were:  648.7,  551.7, 
403.2,  365.4,  199.4,  125.2,  80.5,  45.8,  34.1,  16.4,  6.2,  2.9,  2.4,  1.2, 
and  1.0.  Note  that  there  are  13  unique  eigenvalues  represented  here. 
Utilising  Time-Delay  Steering  on  All  Element  Signals 

If  it  is  beneficial  to  minimize  beam  shifting  of  the  AS  sector 

reference  beam  by  employing  time-delay  steering,  then  it  should  be  even 

more  beneficial  to  extend  this  characteristic  to  the  Interferometer  Beam- 
former  as  well  by  utilizing  time-delay  steering  on  all  element  signals. 
This  does  indeed  prove  to  be  the  case  for  our  SLPF  estimator /processor 
system,  and  we  end  up  with  a  practical  manifestation  of  the  "focussed" 

covariance  matrix  technique  proposed  by  Wang  and  Kaveh  [33.  Essentially, 

the  time-delay  steering  minimizes  equivalent  spatial  bandwidth  spread  to 
such  a  small  value  within  the  AS  sector  that  transversal  filters  are  not 
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really  needed  (assuming  modest  RF  bandwidths).  This  benefit  is  demon¬ 
strated  in  Fig.  14(a)  where  we  plot  typical  adapted  patterns  using  time- 
delay  steering  on  all  element  signals,  with  no  transversal  filters 
present.  All  nine  interferometer  beams  were  used,  thus  giving  us 

A 

a  R  matrix  dimension  of  9,  and  the  total  number  of  snapshots  averaged  was 
N  -  315. 

The  most  remarkable  feature  about  this  technique  is  that  we  are 
getting  an  estimation  performance  fully  equal  to  the  performance 
demonstrated  in  Fig.  13,  but  tith  only  3  DOF  needed  as  compared  to  13  for 
the  TF  configuration.  The  DOF  information  is  contained  in  the 

A 

eigenvalues,  which  computed  for  the  particular  R  of  this  case  as:  400.3, 
345.9,  101.8,  1.4,  1.2,  1.1,  1.0,  0.95,  and  0.91.  Note  that  there  are 
only  3  unique  eigenvalues  represented  here,  which  also  happens  to  be  the 
minimum  number  associated  with  three  spatial  point-sources.  This  clearly 
indicates  a  "focussed”  covariance  matrix  and  would  permit  us  to  reduce  the 
matrix  dimension  down  to  3  DOF,  i.e.,  only  three  assigned  beams  for  this 
particular  case.  Additional  benefits  include  fewer  data  snapshots 
required  and  faster  speed  of  response. 

Time-delay  steering  to  a  narrow  AS  sector  is  a  fundamental  technique 
for  separating  the  spatial  domain  from  the  frequency/time  domain.  When 
combined  with  a  shaped  AS  sector  reference  beam,  it  functions  in  much  the 
same  manner  as  a  narrowband  filter  in  the  frequency  domain,  i.e.,  eitner 
technique  can  separate  the  two  domains  within  the  limitations  of  their 
"passbands"  and  "sidelobe  levels".  To  illustrate  this  point,  the  final 
simulation  example  deliberately  includes  two  additional  20  dB  sources  of 
10J  swept  bandwidth  located  at  —13-5  and  -12.2  degrees  azimuth.  Referring 
to  Fig.  11,  note  that  these  sources  will  suffer  at  least  24  dB  of  attenua¬ 
tion  in  the  reference  beam  sidelobes  and,  therefore,  will  contribute  very 
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little  to  the  cross-oor related  steering  vector  A  of  Eq.  (17).  In 
addition,  they  undergo  an  equivalent  spatial  spreading  beoause  they  are 
about  55  degrees  away  from  the  time-delay  steering  direction  of  42°.  The 
resultant  adapted  patterns  are  shown  in  Fig.  14(b)  where  we  note  only 
secondary  disturbances  due  to  the  two  extra  swept  sources.  If  significant 
correlation  were  present,  there  would  be  a  cancellation  notch  extending 
over  the  indicated  aliased  region. 

An  interesting  aspect  of  this  pt”tioular  case  is  that  the  two  extra 
swept  sources  are  presented  at  full  gain  in  some  of  the  Interferometer 
Beamformer  beams  because  of  the  wide  element  spacing  of  310«i.e.,  we  have 
replication  of  the  set  of  beams  in  Fig.  10  throughout  visible  space.  This 
effect  is  threefolds 

a 

a.  They  contribute  terms  of  high  power  level  in  R. 

b.  They  increase  and  dominate  the  eigenvalues,  which  were  computed 
as:  703.6,  575.8,  419.5,  388.1,  287.6,  192.2,  35.3.  2.3.  and  1.1.  Note 
that  there  are  now  8  unique  eigenvalues  instead  of  3. 

c.  They  have  an  equivalent  spatial  bandwidth  spread  in  the  AS  sector 
because  of  the  aliasing  from  their  actual  location.  This  equivalent 
spread  is  indicated  in  Fig.  14(b). 

The  aliasing  intrusion  into  our  AS  sector  can  be  harmful  if  it  overlaps 
true  sources  in  *  the  sector,  because  then  the  adapted  weights  will  be 
affected  by  the  extra  signal  power  coming  from  "common”  directions,  and  it 
may  be  sufficient  to  distort  or  even  destroy  the  null  in  those  directions. 
VI.  CONCLUSIONS 

We  have  Investigated  the  concept  of  performing  SLPF  source  estima¬ 
tion  in  a  narrow,  steerable,  spatial  sector  by  utilizing  a  large  antenna 
aperture  width  (in  wavelengths),  but  requiring  relatively  few  DOF  from  it 
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via  operation  in  an  interferometer  beamformer  space.  Our  objectives 
included  high-resolution,  the  accommodation  of  moderate  bandwidths,  fast 
response  time,  moderate  cost,  and  ease  of  frequent  calibration.  Computer 
simulations  on  several  multiple-source  examples  have  demonstrated  that  the 
concept  is  indeed  feasible,  provided  true  time-delay  steering  can  be 
employed  on  the  AS  sector  reference  beam.  This  provision  retains  a 
reasonably  stable  spatial  reference  for  each  source  within  the  sector  over 
the  common  bandwidth. 

The  optimum  estimation  system  results  when  time-delay  steering  can 
be  applied  to  all  element  signals,  i.e.,  including  the  interferometer 
elements,  because  this  "focusses"  the  sample  covariance  matrix  for  sources 
which  are  within  the  AS  sector  and  minimizes  the  DOF  required.  Even  the 
transversal  filters  may  be  dispensed  with  under  "focussed"  conditions  if 
the  bandwidth  is  modest.  When  combined  with  a  shaped  AS  sector  reference 
beam,  this  technique  functions  in  much  the  same  manner  as  a  narrow  band 
filter  in  the  frequency  domain,  i.e.,  it  separates  the  spatial  domain  from 
the  frequency /time  domain  within  the  limitations  of  its  "passband"  and 
sidelobe  levels. 
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APERTURE 

ELEMENTS 


Fig,  2  —  Schematic  diagram  of  constrained  beamspace  SLC  system 
with  a  tapped  delay  line  on  each  assigned  beam 
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Fig.  4  —  Typical  adapted  patterns  at  low,  midband,  and  high  frequency  for  a  single  46  dB  point- 
source  located  at  42  degrees,  swept  RF  bandwidth  of  10  percent;  single  SLC  beam  No.  14  with  3  tap 
transversal  filter. 


POWER  IN  DECIBELS 


Fig.  7  —  Typical  adapted  pattern  response  for  eleven  40  dB  sources  distributed  uniformly  over  spatial 
sector  width  of  1.2S  beamwidths  (37.7  to  S0.2  deg.);  three  SLC  beams  Nos.  13,  14,  and  IS  without 
transversal  filters. 


Fig.  8  —  Summary  of  adaptive  cancellation  performance  for  eleven  40  dB  sources  distributed 
uniformly  over  a  spatial  sector;  utilizing  one,  two,  or  three  auxiliary  SLC  beams  without  transversal 
filters. 
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Aperture  Elements 


Schematic  diagram  of  a  SLPF  Estimator  incorporating  a  Sparse-Array  Beamformer, 
a  shaped  Reference  Beam,  and  a  Transversal  Filter  structure. 


SPATIAL  ANGLE  IN  DEGREES 


SPATIAL  ANGLE  IN  DEGREES 

(b)  Five  interferometer  beams  with  a  3-tap  transversal  filter  on  each  beam 

12  —  Typical  adapted  patterns  for  the  case  of  three  16.7  dB  sources  of  10%  bandwidth 
located  at  37,  42,  and  48  deg.,  with  no  equivalent  "overlap”. 
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(b)  AS  sector  reference  beam  steered  via  true  time-delay 

Fig.  13  —  Typical  adapted  patterns  for  the  case  of  three  16.7  dB  sources  of  10%  bandwidth  located 
at  40,  42,  and  44  deg.,  with  significant  spatial  "overlap”.  Five  interferometer  beams  with  a  3-tap 
transversal  filter  on  each  beam. 
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SPATIAL  ANGLE  IN  DEGREES 


(a)  Same  swept- sources  case  as  Fig.  13 


(b)  Swept  sources  case  of  Fig.  13  plus  two  additional  20  dB  sources  of 
10%  bandwidth  located  at  -13.S  and  -12.2  deg.  azimuth 

Fig.  14  —  Typical  adapted  patterns  using  time-delay  steering  on  all  element  signals, 

with  no  transversal  filters  present. 
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APPENDIX  A 


Interferometer  B— former  Weights 

The  beamformer  indicated  in  Fig.  9  may  consist  of  a  Butler  Matrix 
design  C A1 3 »  a  Rotraan  lens  design  [A2],  or  any  other  feasible  implementa¬ 
tion  [A3].  Orthogonality  in  the  family  of  beams  is  desirable  but  not 
essential,  and  low  sidelobe  design  is  not  required.  A  Rotman  lens  would 
offer  the  inherent  advantage  of  incorporating  true  time-delay  beam  forma¬ 
tion  such  that,  in  the  light  of  the  "focussed"  sector  discussion  of 
Section  V,  it  would  probably  be  a  preferred  candidate  in  an  actual  system 
design.  However,  for  the  purposes  of  this  report,  where  it  was  necessary 
to  evaluate  both  phase-steer ing  and  time-delay  steering,  the  Butler  Matrix 
was  a  better  choice. 

A  Butler  Matrix  essentially  sets  up  a  spatial  discrete  Fourier 
transform  in  which  the  transformation  vector  for  a  linear  array  with  equal 
element  spacing  will  have  individual  weights  of  the  form. 


where 


bkm  -(/iHexp 

X  -  (f^)  Sine 
Ao 


(A1  ) 


( A2) 


m  -  beam  index 
k  -  element  index 
K  -  total  number  of  elements 
d  -  element  spacing 
0 

b  -  vernier  beam  adjustment  angle. 

For  our  case  where  K  is  an  odd  number  (K-9)  we  chose  9.  such  that 

b 

X«1 ,  and  this  automatically  adjusts  the  center  beam  to  boresight  position. 
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It  should  be  noted  that  the  element  weights  expressed  In  Eq.  (A1 )  are 
independent  of  the  spaoing  of  the  elements,  except  for  the  vernier  adjust** 
ment  faotor,  X.  The  complete  transformation  matrix  £  then  consists  of  K 
column  vectors  of  K  element  weights,  where  the  kth  oolumn  vector  b^  is 
associated  with  the  kth  bean, 


B 


l  l  I  & 
b21  b22 
b31  lb32 

•  I  • 


•  • 


A 

The  beamformer  output  vector  E  is  expressed, 


(A3) 


(A4) 


where  JE  is  the  column  vector  of  input  signals  from  the  array  elements. 

This  beamformer  is  referred  to  as  an  "interferometer  beamformer" 
because  the  input  elements  from  the  aperture  are  spaced  d-31Q  apart.  The 
implication  of  this  wide  spacing  is  that  the  cluster  of  beamformer  beams 
is  not  unique  and  will  replicate  itself  several  times  throughout  visible 
space.  For  example,  if  we  consider  just  the  center  beam  and  its  boresight 
position,  we  will  get  its  full  uniform  illumination  peak  gain  at  multiples 
of  2it  interelement  phasing, 

2iTY“Sin8.  -  2ni  (A5) 

o 
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or  sinS^  -  1/3  for  1  *  *3,  “2,  -1,  0,  1,  2,  3.  Therefore,  We  center  beam 
replicates  at  the  azimuth  spatial  angles  of  -90°,  -41.8°,  -19.5°,  0°, 
+19.5°,  +41 ,8°,  and  +90°. 

AS  Sector  Reference  Beam 

The  reference  beam  is  intended  to  be  a  spatial  filter  which  passes 
the  AS  sector  but  filters  out  all  sources  that  are  located  outside  it.  In 
attempting  to  approximate  this  ideal  rectangular-shape  filter  without 
getting  into  sophisticated  design  procedures,  the  method  of  separating  two 
sin  x/x  beams  was  utilized.  It  is  readily  shown  that  the  element  phasing 
required  for  separating  the  twobeams  by  A  beamwidths  is, 


,  TfAr  K  +  1  1 

■  rtk  ‘  —!• 


(A6 ) 


and  adding  the  two  beams  results  in  cosine  element  weights  Wk, 


1  r 

Wk  -  ^[e  +  e  ]  -  cos^ 


(A7) 


The  array  parameters  chosen  for  thi3  simple  cosine-type  weighting  were 
K-19  elements  spaced  a  half-wave length  apart,  and  a  beam  separation  of 
A  ■  1 .5  beamwidths  where  a  beamwidth  is  defined  as  the  arc  sin(2/K).  By 
itself,  this  weighting  produces  a  double-hump  pattern  with  a  2.2  dB  dip  in 
the  middle,  so  the  end  element  weights  wer*-  increased  by  +.546  to  flatten 
and  extend  the  central  constant-gain  region  Final  weights  obtained  were: 
-.068,  -.402,  -.164,  +.082,  +.325,  +.546,  +.735,  +.880,  +.970,  and  +1.0, 
where  +1.0  is  the  weight  applied  to  the  center  element.  Figure  11 
illustrates  the  resultant  flat -top  beam  at  midband,  steered  to  a  center 
position  of  42  degrees  azimuth. 
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Figure  A1  illustrates  the  pattern  shifting  that  occurs  at  the  lowest 
and  highest  frequencies  within  a  10J  RF  bandwidth  under  the  two  methods 
for  steering  the  beam,  phase-shifters  and  true  time-delay.  The  phasing 
associated  with  the  above  weights  for  these  two  methods  is  readily  deter¬ 
mined  for  steering  angle  8q  and  half-wavelength  spaced  elements, 


phase  shifters  -  exp(-j  1^) 

(A8) 

time  delays  -  exp(-j(f9/f  )¥.) 

Xf  0  K 

(A9) 

rr  «  « 

where  4>k  -  irsin8o  [k  -  -^-j, 

fQ  is  the  midband  RF  frequency  and  f^  is  the  fcth  frequency  within  the 
bandwidth.  The  essential  difference  is  that  the  phase  shifter  is  a  fixed 
phase  setting,  ¥  whereas  time  delay  offers  a  counteracting  frequency 
compensation  to  the  midband  setting. 

A  final  point  to  be  emphasized  in  regard  to  the  AS  sector  reference 
beam  is  that  its  phase  center  is  displaced  -7.5  \Q  away  from  the  phase 
center  of  the  interferometer  beams.  This  displacement  is  essential  to 
achieving  robust  spatial  estimation  performance  in  our  SLPF  systems, 
because  it  produces  a  phase  sensitivity  to  source  location  in  a  beamspace 
system,  in  addition  to  the  amplitude  sensitivity  derived  from  beamshapes. 
One  consequence  of  this  displacement  is  that  it  must  be  taken  into  account 
in  referring  the  AS  reference  beam  output  signal  to  the  phase  center  of 
the  interferometer  beamformer.  Thus,  it  adds  an  additional  term  to  Eqs. 
(A8)  and  (A9), 
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(A10) 


phase  shifters  -  exp{-j  (¥  k+  a)} 
time  delays  -  exp{-j(f#/f  )[¥  ♦  a]  (All) 

0  K 

where  a  -  2ir(7.5)sin0  . 

o . 
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(a)  Beam  steered  to  42  deg.  azimuth  via  phase  shifters 


(b)  Beam  steered  to  42  deg.  azimuth  via  true  time-delay 

Fig.  A1  —  AS  sector  reference  beam  plotted  at  three  frequencies 
within  10%  RF  bandwidth 
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